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Modelling rapid TeV variability of PKS 2155-304 
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ABSTRACT 

We present theoretical modelling for the very rapid TeV variability of PKS 2155-304 
observed recently by the H.E.S.S. experiment. To explain the light-curve, where at 
least five flaring events were well observed, we assume five independent components of 
a jet that are characterized by slightly different physical parameters. An additional, 
significantly larger component is used to explain the emission of the source at long 
time scales. This component dominates the emission in the X-ray range, whereas the 
other components are dominant in the TeV range. The model used for our simulation 
describes precisely the evolution of the particle energy spectrum inside each component 
and takes into account light travel time effects. We show that a relatively simple 
synchrotron self-Compton scenario may explain this very rapid variability. Moreover, 
we find that absorption of the TeV emission inside the components due to the pair 
creation process is negligible. 
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1 INTRODUCTION 

Among the different models proposed to explain X-ray and 
gamma-ray emission of TeV blazars one disti nguishes two 
classe s: the leptonic and hadronic models (e.g. iKrawczvnskil 
120041 '). This division depends on the assumption one makes 
about the particles that initially carry the energy that is 
to be converted into electromagnetic emission during the 
evolution of the source. The aim of this paper is to explain 
the rapid variability in the simplest possible way in order to 
constrain the physical parameters of the source. Therefore, 
we have decided to use a relatively simple leptonic scenario. 
A very basic leptonic model that can be used to ex- 
plain the high energy emission of TeV blazars assumes 
a spherical source filled with a tangled magnetic field 
and relativistic electrons. This source is thought to travel 
with relativistic velocity at a distance of less than 1 
pc from the central black hole and is usually assumed 
to be homogeneous. The electrons inside the source pro- 
duce synchrotron emission and also up-scatter part of 
this emission to TeV energies. This is the well known 
synchrotron self-Compton process (SSC) that is often 
applied to desc ribe X-ray and gamma- r ay observations 
of bl a zars (e.g. [ Blo om & Marschcr 1996, Ghisellini ct al. 
19961 . Ilnoue fc Takaharal Il996l , iM a stichiadi s fc Kirkl Il997l. 



rectly from observations (e.g. Bednarek fc Protheroel Il997l . 



Krawczvnski et al.l 120 00). The simplicity of this model al- 



lows to constrain some important physical parameters of 
the source (e.g. magnetic field strength, Doppler factor) di- 

* E-mail:kat@astro. uni.torun.pl 



iTavecchio et al.lll998l . lKatarzvnski et al.ll200ll ). On the other 
hand, some uncertainty remains in such estimations because 
the model is too simple; it does for example not take into 
account significant processes such as the particle evolution 
inside the source or light crossing time effects (hereafter 
LCTE). Therefore it is very important to analyse not only 
the emission observed at a given time but also the evolution 
of the emission, especially during the periods of activity. 
This may provide additional constraints for the models and 
result in a more precise estimation of the physical parame- 
ters. 

Several leptonic models have been proposed to ex- 
plain the time-dep e ndent emission o f TeV blazars (e.g. 
iDermer et~aH Il997l. iKirk et all Il998t ICoppi fc Aharonianl 
19991, IChiaberge fc Ghisellinil 1 19991, iKusunose et al.l teo'obT 
iKataoka et al.l l2000. So kolov et aT1l2004h . So far, no model is 
able to precisely describe SSC emission of an inhomogeneous 
or even a homogeneous source. The main problem appears 
in the description of the synchrotron radiation field inside 
the source that is inhomogeneous even in a sph erically ho- 
mogeneous slowly evolving source l|Gouldlll979l ). Moreover, 
this radiation field depends on the particle energy distribu- 
tion and the particle energy distribution depends in turn on 
the radiation field. This is certainly only true if the energy 
density of the radiation field is comparable to or larger than 
the energy density of the magnetic field, but this condition is 
fulfilled in TeV blazars. Finally, very rapid variability time 
scales are usually assumed to be comparable to the light 
crossing time at the source. This means that LCTE must be 
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taken into account in calculations of the inverse-Compton 
scattering. In other words, the radiation field at a given time 
and a given position inside the source depends not only on 
the local physical conditions, but also contains photons from 
other parts of the source. These photons were already cre- 
ated at an earlier time and have propagated to the observed 
position. We will call this process internal LCTE. This effect 
is especially important for the calculation of the inverse- 
Compton emission and somewhat less for the synchrotron 
emission (although the radiation field has generally also an 
impact on the evolution of the electron energy spectrum). 
The internal LCTE was tak en into accoun t for th e first time 
in the model proposed by ISokolov et al.l ((200J). However, 
this model assumes dominance of the energy density of the 
magnetic field over that of the radiation field. In this case 
the impact of the radiation field on the electron energy spec- 
trum is negligible. 

A second LCTE, which we call external, is equally im- 
portant for both the synchrotron and the inverse-Compton 
emission. The external observer receives at a given time the 
emission produced by different parts of the source at dif- 
ferent times in the comoving frame of the source. This is 
due to the different travel times that photons from different 
parts of the source need to reach the observer. If the emis- 
sion level of the source does not change, the observer will 
always receive the same amount of radiation; if it changes 
on long time scales, differences in photon travel times can 
be neglected. However, if there is a change of the emission 
level with a duration shorter than the light crossing time of 
the source, the different travel times of the photons must 
be taken into account. On e of the first attempts to ac count 
for this effect was made bv lChiaberge fc Ghisellinil (|l999T ). It 
should be noted that in their model they neglected internal 
LCTE, which is a good enough approximation under certain 
conditions, as we will explain below in the description of the 
model. 

Another important problem in time dependent SSC 
modelling is the description of the particle acceleration and 
evolution of the particle energy spectrum. It is widely be- 
lieved that the particles are accelerated by shock waves cre- 
ated by colliding components of a jet. In such collisions, 
some fraction of the kinetic energy related to the bulk mo- 
tion of the components is transformed into random kinetic 
energy of the particles. This energy is then radiated away. 
The acceleration and evolution of the spectrum are usually 
described by the k inetic equation, w hich is a partial differen- 
tial equation (e.g. iKardashev Ill962f ). There are two general 
approaches to describe the acceleration in the equation. 

In the first approach, the acceleration term in the equa- 
tion describes the process inside a shock wave. T his ap- 
proach was used for example by iKirk et al.l (|1998l ) , where 
the evolution of the source was represented by two kinetic 
equations. The first equation described the acceleration in- 
side the shock and the second the evolution of the particle 
spectrum in the downstream region of the shock. Particles 
were escaping into this region after the acceleration, cre- 
ating a source for electromagnetic emission. However, this 
is a rather complex approach and was applied only for the 
synchrotron emission. 

The second approach is simpler. There is no acceleration 
term, i.e. no description of the acceleration process. Instead, 
an "injection term" defines directly the result of the accel- 



eration. One assumes a particle energy spectrum created by 
a hypothetical shock and describes it with a few free pa- 
rameters. T his approach was used fo r exam ple in the model 
proposed bv lChiaberge fc Ghisellinil |l999j), where different 
types of injection were investigated. 

Since there is no perfect solution for modelling of vari- 
ability in the framework of the leptonic SSC scenario, we 
have chosen here the relative ly simple model proposed by 
iChiaberge fc Ghisellinil (Il999l ). The model, as we will show, 
provides very reasonable results while using only a minimum 
number of free parameters. 



2 RAPID VARIABILITY 

High energy variability of TeV blazars has been observed 
many times by orbiting X-ray experiments, as well as ground 
based gamma-ray telescopes. One of the most spectacular 
states of high activity was observed in Mrk 501 in 1997. The 
emission level and the spectrum of the emission changed dra- 
matically in a p e riod of a few days dCatanese et al.lll 997. 
Pian et al.l Il998l . iDiannati et al.l Il999l . iKrawczvnski et al.l 



20001 ). Very rapid variability at time scales of only a few 



hours was first observed simultaneously in X-rays and TeV 
gamm a rays in Mrk 421 and reported by iMaraschi et al.l 
(1999). Today Air Cherenkov Telescopes like H.E.S.S. or 
MAGIC achieve major breakthroughs with minute-scale ob- 
servations in the VHE range. 

PKS 2155-304 has been intensively observed during the 
last years by different instruments, showing variability at 
time s cales from several m onths to a few hours. For ex- 
ample iGiommi et al.l (|l99a ) observed X-ray activity of this 
source with a variability time scale of a few hours. Simi- 
lar time scales o f the X -ray activity were also r eported by 
IChiappetti et all 119991) andlEdelson et al.1 (J2001). Observa- 
tions made by iTanihata et al.l (|200lh show X-ray flux vari- 
ations at time scales of about one day. In the TeV range, 
the sou rce has been detected by the Durham Mk VI tele- 
scopes l|Chadwick et al.lll999T ). This detection remained un- 
confirmed until observatio ns were made by the H .E.S.S. ex- 
periment in 2002 and 2003 (lAharonian et alfcOOSl ). Variabil- 
ity on daily time scales was observed during that period. 

Recently (28th of July 2006, i.e. MJD 53944), a huge 
TeV outburst of PK S 2155-304 has been observed by the 
H.E.S.S. experiment (|Aharonian et al.ll2007h . The TeV emis- 
sion level of the source changed by a factor of about 20 
within the observational period of about 1.5 h. Variabil- 
ity with a time scale of only a few minutes could be ob- 
served, which presents the fastest variability ever seen in 
a blazar. Unfortunately, there are no simultaneous X-ray 
observations that could give us strong constraints on the 
emission process during this first huge outburst. After this 
outburst, the source was being monitored by several exper- 
iments and especially almost continuously by the UVOT, 
XRT and BAT exp eriments on board the Swift satellite 
l|Foschini et al.l I2007T ) for a duration of several days (from 
29th of July to 22nd of August 2006, i.e. MJD 53945 to 
MJD 53969). The comparison of the Swift d ata with the 
observations made by H.E.S.S. at that time l|Rauel 120061 ) 
shows a rather mi nor change of the U V^X-ray emission 
level (see Fig. 3 in iFoschini et al.l I2007T ) and a huge varia- 
tion of the TeV emission. Compared to most observations 
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of blazars, this is a rather unexpected behaviour; however, 
even more extreme events, so-called orph an TeV flares, have 
already been observed in other sources (JKrawczvnskill2004 . 
iBlazeiowski et a l, 2005). We will in the following try to ex- 
plain this rapid TeV variability event in order to constrain 
the physical parameters of the source. 



3 THE MODEL 



Here we use the model proposed by IChiaberge fc Ghisellinil 
1999. They assume that the source of VHE emission is cre- 
ated by a shock wave that accelerates the electrons. The 
shock is perpendicular to the jet symmetry axis and moves 
with constant relativistic velocity (v' s ). The particles accel- 
erated by the shock up to extremely relativistic energies es- 
cape into the downstream region of the shock, where they 
generate most of the synchrotron and IC radiation. In other 
words the shock is injecting relativistic particles into some 
volume where most of the particle energy is dissipated by 
radiative processes. Therefore, this volume is equivalent to 
the source volume. For sake of simplicity, the cross section 
of the shock front perpendicular to the jet symmetry axis 
is assumed to be square, of scale R. The thickness of the 
shock front (parallel to the jet axis) is r <C R- The shock is 
accelerating the particles and thus creating the source for a 
duration of t CI = R/c. Such an assumption eliminates the 
duration of the acceleration process as a free parameter and 
yields a final volume of the source of R 3 . The evolution of 
the source is also simulated after the injection phase for a 
few crossing times up to almost complete decay of the X-ray 
and the 7-ray emission. 

The source volume is divided into 10x10x10 smaller 
cells and calculations are performed in time steps of At' = 
O.lR/c. This allows to take into account the external LCTE, 
which is the main advantage of the model. In the first step, 
at time t' = At' , only a narrow region just after the shock 
(at the distance from to v' s At') is filled up by particles. In 
the next step (t' — 2At') the particles injected during the 
first phase are somewhat "older" which means they have 
already lost some fraction of their energy. They are located 
at a distance between v' g At' and 2v' s At' . It should be noted 
that, when solving the kinetic equation that describes the 
evolution of the particle energy, we assume that the injection 
process works only during a time ti n j = O.lR/c, necessary 
to create a cell. Moreover injection is only active in the cells 
that are next to the shock front; for other cells we calculate 
only radiative cooling. The shock is creating continuously 
new cells filled by "fresh" particles, whereas old cells are 
moving systematically to larger distances. However, what 
is observed in the comoving frame is significantly different 
from this due to LCTE. The observed spectrum is produced 
by the electron distribution at different stages of evolution. 
Therefore, to obtain the total emission, it is necessary to 
sum up the different contributions of different cells in an 
specific way, as described in details in the original paper. 

The source is observed at an angle 9 — 1/F [rad] , where 
F is the Lorentz factor that describes bulk motion of the 
source. This means that in the source's comoving frame the 
emission is observed at 90 degrees with respect to the jet 
symmetry axis and the shock velocity vector. For such an 
angle the external LCTE has the strongest impact on the 



evolution of the observed emission. Note that for zero degree 
the external LCTE is negligible. 

IChiaberge fc Ghisellinil (|l999l ) tested a few different 
types of injection spectra, here we use only the very simple 
single power-law injection Q("j) = Qo7 _ " for 7mm ^ 7 ^ 
7max, where 7 is the Lorentz factor that describes the en- 
ergy of the particle. We neglect also a possible escape of the 
particles from the source. In principle the particles may es- 
cape into regions of the jet where the magnetic field strength 
is significantly lower in comparison to the field inside the 
source. Therefore, synchrotron emission from such regions 
should also be significantly lower and thus negligible. How- 
ever, if the radiation field outside the source is high enough, 
the escaped particles may still produce a significant amount 
of emission through IC scattering. This is a quite complex 
scenario with additional free parameters not yet constrained 
by the data. Therefore, for sake of simplicity, we neglect here 
any possible escape. Morover, we have also introduced some 
small improvements to the model. We apply a more precise 
description of the IC emissivity , which uses the Compton 
kernel computed bv lJonesl l|196g) . The kernel describes scat- 
tering on an isotropic distribution of soft photons, consid- 
ering the full Klei n-Nishina cross section in the head-on ap- 
proximation (e.g. Inoue fc Takaharal 1 19961 . ISauge fc Henri I 
|2004 iModerski et all 120051 ). The cooling rate (7) due to 
IC scattering is calculated more precisely than in the origi- 
nal model, using the Compton kernel mentio ned above (e.g. 
ISauge fc He"nriTl2004 fModerski et all l2005h . We calculate 
the absorption of the TeV emission inside the source due to 
electron-po sitron pai r prod uction using the approximatio n 
derived by ISvenssor] l| 19871 ) (e.g. Ilnoue fc Takaharal I1996J ). 
The absorption is calculated separately inside each cell and 
absorption caused by the surrounding cells during the prop- 
agation of the emission from a given cell toward the observer 
is also taken into account. To correct the intrinsic TeV spec- 
tra for the IR e xtragalactic absorpti on, we use an optical 
depth deri ved by Stecker et al.l (|200q ) with the appropriate 
correction IjStecker et al.l 120071 ). which was calculated in a 
ACDM universe for h = 0.7, Q A = 0.7 and Q. m = 0.3. 

There are several clear advantages to this model. The 
most significant one lies in the reduction of the number 
of free parameters needed to fully describe a single flar- 
ing event. For single power-law injection we need in prin- 
ciple four free parameters: minimum and maximum energy 
(7min,7max), density of the injected particles (Qo) and in- 
dex of the power-law (n). However, precise values of 7 m i n 
and 7 ma x are not crucial for the model. To provide a sig- 
nificant amount of low energy synchrotron photons for the 
IC scattering, it is indeed sufficient to keep 7mm relatively 
small (e.g. 1 ^ 7mm ^ 10 3 ). On the other hand, the value of 
7max inside each cell at a given time is precisely calculated 
according to the cooling conditions inside the source. There- 
fore, it is sufficient to assume 7 max high enough to be able 
to produce TeV gamma-rays. Simulating the emission of six 
different sources we then use 7 m in = 1 and 7max = 10 6 . The 
index of the injected spectrum, n = 2, is the same in all 
our calculations. Three more free parameters are required 
to completely describe a single source: the source extension 
(R), the magnetic field strength inside the source (B), and 
the Doppler factor (<5). This results in 7 free parameters, a 
number identical to the one required in the simple station- 
ary scenario (see the introduction). However, if this simple 
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scenario uses a broken power-law particle energy distribu- 
tion (this is required in most cases) then the number of 
free parameters increases to 9. Moreover, it is not clear how 
such a broken power-law energy distribution could be cre- 
ated when for example the spectral index above the break 
is > n + 1. On the other hand, the model we use requires 
only a single power-law injection to provide a self-consistent 
explanation for the observed X-ray spectra that one could 
approximate by broken power-law functions. Such broken 
synchrotron spectra are created by external LCTE. This is 
another advantage of the model and this shows also that 
the observed index of the synchrotron emission may have 
nothing to do with the index of the particle energy distribu- 
tion inside the source. Since 7 m i n , 7 max and n have the same 
values in all our calculations and minimal impact for the 
result of the modelling, we are left with only four important 
parameters (Qo, B, R, S) that describe the evolution of a sin- 
gle source. In comparison to the nine parameters required 
by the simple stationary scenario, this is a quite simple de- 
scription for a single activity event. This is important when 
one has to explain several flares one by one. 

The model does not describe internal LCTE, which is 
the main drawback of this scenario. However, this effect is 
not very important in the case where single power-law par- 
ticle energy distribution is injected. After the injection, the 
particle energy spectrum inside the cells is modified only by 
radiative cooling that is systematically reducing 7 max . This 
means that the synchrotron emission produced by the most 
energetic particles is also systematically reduced. However, 
this emission contributes little to the IC scattering due to 
the Klein-Nishina restriction. On the other hand, the syn- 
chrotron emission produced by the particles with medium 
and low energies remains unaffected by the cooling for a rel- 
atively long time. Thus the synchrotron emission produced 
by those particles remains equal and constant inside most 
of the cells and is largely responsible for the IC scattering. 



4 RESULTS 

There are no X-ray observations simultaneous with the rapid 
variability observed by H.E.S.S. on the 28th of July 2006. 
Therefore, to constrain the model we decided to use the data 
obtained by Swift during the several day long campaign of 
observations that was triggered by this rapid TeV activity. 
The X-ray light curve and spectrum does not show a signifi- 
cant change of the emission level, whereas TeV observations 
(Raue 2006) shows a dramatic change of the TeV emission 
level during this campaign. Thus, we decided to simulate an 
event where X-ray emission remains almost constant during 
the TeV flaring activity. As we already mentioned such ac- 
tivity was already observed at least two times in other TeV 
sources. 

However, such a situation is problematic for SSC 
modelling, where the synchrotron emission that provides 
low energy photons for the IC scattering should usually 
vary as much as the IC emission. Only in the quite un- 
realistic case where only the particle density inside the 
source is varying may the IC emission increase or decrease 
as the square of the synchrotron emission (Fic(v,t) oc 
f-Fsvnch(^i^)1 )■ However, this anyw ay may be diluted by 
the LCTE (JKatarzynski et al.l [20051 '). The simple solution 



Table 1. The most important parameters used for the modeling, 
where c\ indicates background component of the jet and C2^6 
indicate parameters that describe foreground components. The 
other parameters: n = 2, 7 m i n = 1 and 7 max = 10 6 are the same 
for each component. Note that we specify particle density inside 
source (ti n jQo) instead of the injection rate (Qo) to give possibil- 
ity for easy comparison between our calculations and results that 
may give simple stationary SSC scenario. 



name 


Cl 


C2 


C3 


C4 


C5 


C6 


unit 


S 


20 


30 


30 


30 


30 


30 




H 


3 10 3 


5.2 


3.5 


4 


2.3 


6 


10 14 cm 


B 


0.02 


0.1 


0.1 


0.1 


0.1 


0.1 


G 


*injQ0 


8 10" 6 


3.2 


7.95 


6.8 


12 


2.45 


10 7 cm -3 



for this problem is to assume that the X-ray emission is 
dominated by a relatively large component from a jet that 
provides almost constant emission on a time scale of several 
days. Since the density of such a large component is rela- 
tively low, the level of the TeV emission must also be neg- 
ligibly low. The existence of such a component is strongly 
supported by the observations made by Swift, where for sev- 
eral days a nearly constant emission level was observed. The 
H.E.S.S. observations show that, in addition, relatively small 
components can appear inside the jet. These components 
may not be strong enough to be dominant in X-rays, but 
may be very compact and therefore very strong in the TeV 
range. 

We simulate in our model one large jet component (we 
call it background component because it is not visible in the 
TeV range), and five small components (here called fore- 
ground components) dominant at TeV energies. The main 
difference between the background component and the fore- 
ground components appears in the estimate of the size and 
the particle density. Moreover, the small components are 
likely located at a distance of less than 1 pc from the cen- 
ter, whereas the large component should be located further 
downstream in the jet, a few parsec from the center. 

In Fig. I we show the light curve and average spectrum 
obtained by H.E.S.S. on MJD 53944 as well as the spec- 
tra obtained by Swift a few days later. We also show in 
this figure the results of our calculations, where rapid TeV 
variability is well reproduced by the superposition of the 
IC emission of the five foreground jet components (denoted 
C2->e)- Note that the IC emission of the background com- 
ponent (ci) is negligible in the TeV energy range. On the 
other hand, this component dominates the emission in the 
X-ray range. We show a theoretical light curve calculated 
in the Swift energy range 0.3-10 keV, where emission from 
the background component appears to be almost constant 
in the relatively short time interval (2.5 h) presented in the 
figure. The total duration of the background component ac- 
tivity, which is calculated in the same way as the activity 
of the foreground components, is about 15 days. This long 
term activity starts 10 days earlier than the first foreground 
component (c.2) flare. Finally, we show synchrotron and IC 
spectra calculated for two arbitrarily selected times (MJD 
53944.035 and 53944.049) that correspond to medium and 
high TeV emission levels. Note that the calculated spectra 



Modelling rapid TeV variability of PKS 2155-304 5 



3.5*10 -9 - 



-9 



3.0*10' 



eg 2.5*10 -9 - 



1*10 -10 



U. d 0*10° 



I I I I I I I I I I I I I I I 

c) 

H.E.S.S. MJD 53944 ■ 
Swift 
MJD 53946 and 53949 



1 c, + c, 




c, + c, 



0'' 'V 0-' 



-8 IT- 1 -8 




I I I I I I I I I I I I I I I 
d ) H.E.S.S. MJD 53944" 



Swift 

MJD 53946 and 53949 




0' y\ w 

i A i I i I A I i^ 




12 14 16 18 20 22 24 26 28 

log.Cv) [Hz] 



H.E.S.S. E>200 GeV 




b) 







0.0 




1.0 1.5 

time from MJD 53944.0 in hours 



2.0 



2.5 



Figure 1. Activity of PKS 2155-304 observed by H.E.S.S. llAharonian et al.ll2007h & Swift JFoschini et alj |2007l in July and August 
2006 and results of our modelling. Panel a) shows the rapid variability observed by H.E.S.S. and the calculated light curve (thick solid 
line), where the contributions from simulated jet components (c2— >6) a re indicated by thin solid lines. A theoretical light curve has also 
been calculated in the X-ray range (panel b). Here, an extended, slowly evolving jet component (ci) dominates the emission over the 
other jet components. Sub-panels c) and d) show the observed average spectrum in the TeV range and two spectra obtained by Swift 
a few days later. In the sub-panels, we also show spectra calculated for two arbitrarily chosen times (MJD 53944.035 and 53944.049) 
corresponding to medium (c) and high (d) TeV emission levels. The solid lines in the sub-panels show the total emission; the dashed 
lines indicate contributions from individual components. One can see that the background component (ci) dominates in the X-ray range, 
whereas the foreground components (here C2 and C4) dominate the emission in the TeV range. Calculated not absorbed intrinsic TeV 
spectra are indicated by doted lines. 



are not well constrained, since we have only an average ob- 
served TeV spectrum and no simultaneous X-ray observa- 
tions. Detailed values of the physical parameters used for 
the calculations are given in Tab. 1. Note that in order to 
precisely explain the TeV light curve we had to modify only 
the sizes of the small sources that is constrained by the ob- 
served variability time scales and the particle injection rate 
to explain correctly the level of the emission. 

Rapid variability requires a relatively small and dense 
source, which may be optically thick for TeV emission due 
to pair production inside the source. The simplest solution 
for this problem is to assume a relatively large value of the 
source Doppler factor (5 > 50), which decreases the ob- 
served variability time scale and increases the observed emis- 



sion level. However, such a large value of the Doppler factor 
might also significantly increase the scattering of electrons 
inside the source on the external ambient radiation field. 
Therefore, the external IC scattering might become the dom- 
inant process that produces most of the TeV emission. Suc h 
a scenario was recently proposed bv lBegelman et al.l IJ2007I) . 
However, the efficiency of the external IC scattering depends 
on the intensity of the radiation field and also on the dis- 
tance from the center where the source is created. It may be 
possible that the intensity is too low and the distance is too 
large for this process to be important, as has usually been 
assumed in the modeling of the TeV blazar activity up to 
now. 

An alternate scenario to the external IC scattering was 
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proposed recently by iGhisellini fc Tavecchiol i|2008r h They 
assume that a fast jet (S < 30) contains even faster mov- 
ing compact sources (5 > 50). The radiation field produced 
by the jet is amplified in the comovig frame of the compact 
sources due to the difference in velocities, which increases 
the IC emission produced by these sources. However, mea- 
surements of the motion of parsec-scale jet components in 
PKS 2155-304 at radio frequencies suggest a v alue of the 
Doppler factor of a few (|Piner fe Ed wards! [2004 ) . 

We have carefully calculated the absorption due to pair 
production inside the source using a moderate value of the 
Doppler factor (5 = 30). We found this process to be negligi- 
ble in our approach. There are two reasons for this. First, the 
variability time scale that we are trying to explain is still rel- 
atively long (10-15 minutes), and second, more importantly, 
in all our calculations of the emission from small sources, the 
synchrotron emission level is significantly lower than the IC 
radiation level. It should be noted that we have also verified 
that by trying to explain a variability time scale of about 5 
minutes, as well as by trying to keep vF syn , x-ray — ^fic,TeV 
that is frequently observed in TeV blazars, we may obtain 
significant absorption of the TeV emission. Our results show 
that the classical SSC approach used frequently to explain 
emission of TeV blazars can also explain rapid TeV variabil- 
ity if we assume that the X-ray emission is dominated by an 
extended, slowly evolving source. 



5 CONCLUSIONS 

We have shown in this paper that the rapid TeV variability 
of PKS 2155-304 can be well explained using a standard SSC 
approach while taking into account the particle evolution 
and the external LCTE. In our approach the internal LCTE 
is of minor importance. The model we use requires in prin- 
ciple four free parameters to describe a single activity event. 
Therefore, we may compare the physical parameters derived 
from such modelling with the parameters d erived form pre- 
vious estimations of blazar emission (e.g. iTavecchio et al.l 
Il99a ). This shows that our parameters are similar to other 
results (e.g., B ^ 0.1 G and 8 — 20 — > 30); only the size 
(R ~ 5 x 10 14 cm) of the small sources is about one order 
of magnitude smaller than the values estimated so far. 
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